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Abstract. Generic quantum states in the Hilbert space of a many body system 
are nearly maximally entangled whereas low energy physical states are not; the 
so-called area laws for quantum entanglement are widespread. In this paper 
we introduce the novel concept of entanglement susceptibility by expanding the 
2-Renyi entropy in the boundary couplings. We show how this concept leads 
to the emergence of area laws for bi-partite quantum entanglement in systems 
ruled by local gapped Hamiltonians. Entanglement susceptibility also captures 
quantitatively which violations one should expect when the system becomes 
gapless. We also discuss an exact series expansion of the 2-Renyi entanglement 
entropy in terms of connected correlation functions of a boundary term. This is 
obtained by identifying Renyi entropy with ground state fidelity in a doubled and 
twisted theory. 
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1. Introduction 

The quantum-mechanical state space of a many-body system is huge, but an 
exponentially large fraction of it is inhabited by fairly unphysical states [Tj. Roughly 
speaking one may say that, since we live in a relatively cold universe ruled by local 
quantum held theory, the vast majority of states one encounters in physically relevant 
situations are highly non generic as they are low-energy states of local Hamiltonians. 

For example it has been known for some time now that generic random quantum 
states are nearly maximally entangled P] and that this generic property implies a 
so-called volume law for the entanglement entropy S. More specifically partitioning 
a system into two regions A and B with Na and Nb spins respectively and a with 
a common boundary dA, the entanglement between A and B is 0(min{|^4|, |£>|}). 
At variance with this generic prediction, area laws i.e., S — 0(\dA\). for quantum 
entanglement (as well as other information theoretic quantities [3]) are ubiquitous 
in physics: from black hole thermodynamics to quantum entanglement in zero- 
temperature many body systems (see [I] for a review). 

Moreover also violations of these laws have a fundamental physical meaning, as 
they signal the divergence of an underlying correlation length e.g., quantum criticality 
[5] or give rise to sub-leading terms that allow one to detect hidden topological order 
[6] . Ensembles of physical random states have recently been shown to feature an area 
law in the typical case [7J. Finally, the very fact that area laws hold for many systems 
of interest i.e., the relevant states are just slightly entangled, is the key feature at the 
grounds of efficient classical simulations algorithms for quantum systems [8]. 



While locality and finite range correlations for gapped system certainly provide an 
intuitive ground for understanding the origin of area laws, proving them in a rigorous 
fashion is often quite hard [9] . For the sake of concreteness let us focus on the ground 
state of a local Hamiltonian: what makes it different from a generic quantum states? 
More specifically, why is an area law for entanglement obeyed as opposed to a volume 
law? 

The goal of this paper is to show how one can produce a surprisingly simple 
argument for the emergence of bi-partite entanglement area laws for low-energy states 
e.g., ground states of gapped systems with local Hamiltonians. The argument we 
are going to discuss is both quantitative and general and is based on elementary 
perturbative expansion of the 2-Renyi entropy with respect the boundary Hamiltonian 
terms. This allows one to introduce a natural notion of entanglement susceptibility 
that in turn, and in spite of its simplicity, is able to unveil scaling behavior. 
Remarkably entanglement susceptibility also shows in which ways area laws might 
be violated for gapless systems and excited energy eigenstates. 

In this paper we also discusses an intimate connection between bi-partite 
entanglement and ground state fidelity [10] as well as a relation between entanglement 
susceptibility and fidelity susceptibility (TTJ [12] . Fidelity and entanglement are two 
of the most fundamental quantum information tools and both can be exploited 
to study quantum phase transitions |13] . Entanglement shows that at criticality 
unique quantum correlations become long-ranged and fidelity witnesses the ensuing 
orthogonality catastrophe in the Hilbert space. We find it conceptually rewarding that 
these two phenomena appear now as two sides of the same coin. 
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2. Entanglement susceptibility 

Let us consider a spin system with state-space Ha := <&i£\hi = (C d )® ' A ' and let 
(A, -B) be a bi-partition of the set of vertices A such that, with obvious notation, 
H\ = Ha ® Hb- Any local Hamiltonian H = J2xcA^ x over ^ A ^ x e B{Hx) 
can be written as H = Ha + Hb + Hg where H A /b = ExcA/b ^ and Hg = 
ExnAns#0^ x - We are now interested in the eigen-state properties of the following 
family of Hamiltonians associated to H: 

H(X) = H A +H B +XH d , (AG [0,1]) (1) 

Of course H(0) := Hq describes the situation in which the two regions A and B 
are decoupled whereas H(l) is nothing but the original Hamiltonian where boundary 
coupling Hg is assumed to be non-trivial. The eigenstates of Hq have the factorized 
form \pa)\pb) with eigenvalues E PA + E PB , the non-degenerate ground state being 
l^o) = IO^IOb) with energy Eq. Let |^o(A)) denote the ground state of (fTJ) and pa '■= 
TrBl'I'o(A)) ('I'o(A)l the reduced density matrix of the region A. In order to quantify 
the entanglement between A and B we use the Renyi 2-entropy S2 — ~ logTrp^. The 
first result we would like to present is: 

S 2 (X) = 2X 2 XE + 0(X 3 ) 

(paPb\H9\0a0b) " 



Xe ■■= Yl 

Pa,Pb>1 



E PA + E pB — Eq 



(2) 



In the remainder of the paper we will refer to this quantity as entanglement 
susceptibility [14J. The derivation of ([2]) is straightforward based on Rayleigh- 
Schrodinger perturbation theory and is given in the Appendix. Now a couple of 
simple qualitative comments: 

1) The entanglement susceptibility @ depends explicitly on the Hamiltonian term 
Hg. This dependence makes manifest the intuitive fact that entanglement between A 
and B is caused by interactions terms supported on the boundary between them. 

2) The structure of ([2]) clearly indicates that large gaps (density of states) above 
the ground state energy may result in small (large) value for \e- It should also be 
clear, that vanishing gaps may induce singular behavior of xe depending on the 
behavior of the corresponding matrix elements of Hg. This last fact establishes a direct 
and transparent connection between the behavior of entanglement susceptibility and 
quantum criticality. 



3. The area law bound 



To obtain the area law from (J5J we introduce the so-called fidelity susceptibility 
HUE] Xf ■■= -d 2 logJ"(A)<9A 2 | A=0 , where F(X) := |(* (A)|* )| is the ground state 
fidelity. T is the main ingredient of the so called fidelity approach to quantum critical 
phenomena [10] and has been thoroughly investigated over the last few years. The 
basic fact is that: 

XE < XF (3) 
Indeed xf , for the Hamiltonian family Q at A = 0, is given by 
T, PA+PB >i I {P e PB +e o1 °-e B ) I 2 M that has the ver y same structure of the term we are 
willing to bound but without the restriction that both pa and ps are larger than zero 
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|16j . The bound Q is not always tight. In fact, in order to have potential singularities 
in xe the gaps of subsystem A and B have to vanish. On the contrary for \f the 
closure of the gap in any of the two subsystem e.g., say the largest one is approaching 
criticality, may be enough to generate a divergent behavior. Now we show how the 
bounds proven in [T2] can be adapted to the current case to provide an area law. One 
can write \f = || Goi?a l^o) || 2 where Go = J2 P >o \P){p\/(Ep " ^o), and observe that 
for gapped systems A := E\ — E > and ||Go|| : = sup||^|| =1 ||GoV>|| < A -1 from 
which one has 

Xf < ||G || 2 ||Q# 9 |*o>|| 2 < A- 2 (* |iW/ 9 |*o} 

= A" 2 «* |fl||* ) - (*o|^9|*o> 2 ) • (4) 

This last factor is nothing but the connected correlation function for the boundary 
operator Hg in the ground state ^Jq. Thanks to the exponential decay of correlations 
for gapped ground states of local Hamiltonians [17] the connected correlation above 
scales as the support of the involved operators. More precisely if Hg = J^jedA hj then 
\(HqHq) c \ < £ i d ^ 1 max j £ g a \\hj\\ 2 \dA\, where £ is of the order of the finite correlation 
length of the system measured in units of some microscopic length scale e.g., lattice 
spacing. Therefore, using ([3]), one gets 

X^^maxllM 2 ^- 1 !^!, (5) 

In words: the entanglement susceptibility for a gapped (non degenerate) ground state 
of a local Hamiltonian obeys an area law. We would like to emphasize that this 
argument shows that violations of the area law i.e., of the bound ([S]), are a sufficient 
condition for criticality (for local Hamiltonians). On the other in [TJ] we showed that 
in the case in which Hg has a large scaling dimension criticality does not necessarily 
imply super-extensive scaling of fidelity susceptibility. In the present case this suggests 
that one may have gapless systems that still obey an area law. Such phenomena have 
already been pointed out in the context of projected entangled pair states in [18] or for 
topological insulators which have gapless edge modes but a finite correlation length in 
the bulk PI]. 

Another interesting point is the validity of these results for excited Hamiltonian 
eigenstates. The derivation of @ shows that as long as the excited state is energetically 
well separated from the rest of the spectrum ([2]) holds true. Moreover if the eigenstate 
is also clustering i.e., connected correlations functions of local operators decays 
sufficiently fast, also the derivation of §5§ is valid and hence an area law fulfilled. 
The way this may fail is a high (unperturbed) density of states that makes the excited 
state (quasi-)degenerate and/or a lack of clustering. These remarks are consistent 
with the findings of [2D]. 

4. Example: Quasi-free systems 

The area law prediction of Eq ([5]) for gapped local Hamiltonian can be confirmed by 
explicit, and independent, calculations for exactly solvable models. In this section, for 
the sake of illustration, we discuss two cases: 

a) Free fermions on a graph A with quadratic Hamiltonian H = ^ i jeA Zijc\cj. 

The a's are canonical fermionic operators i.e. {ci,Cj} = &y, {ci,Cj} — and 
Z = (Zi,j)ij is a real symmetric |A| x |A| matrix, local on A [27]. The analog of 
equation ([T]) in this case is given by Z{\) = Z + X5Z where SZ is with non vanishing 
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entries on a 0(|cM| x \dA\) sub-block (rank(<5Z) = 0(|cL4|). The gapfulness assumption 
in this single particle picture amounts to the invertibility of Z i.e., the minimum Z 
singular value A is bounded away from zero. The fidelity susceptibility is known to 
have the form \f = H^V^H! where T(A) := Z\Z\~ X is the unitary part of the polar 
decomposition of Z and ||X|| 2 := VTrXtX [29J. We can now use the bound (see 
chapter VII.5 of [28]) 

\\T A - T B \L < r \\A - B\L 

2 ~ WA-^- 1 + \\B- 1 \\- 1 2 

(with ||X|| indicating the operator norm of X), for A, B invertible operators with 
unitary factors Ta,b in their polar decomposition. Setting A = Z and B = Z + SZ in 
the above inequality we get ||5T|| 2 < 2/(A + A') ||<5Z|| 2 with A' indicating the smallest 
singular value of Z + SZ. Moreover \\SZ\\j < r&nk(5Z)\\5Z\\ 2 with \\SZ\\ = 0(1). 
Thanks to ^ this allows one to conclude 

Xe < (A+ 4 A0 2 rank (^)II^H 2 = ^OQdA\). (6) 

b) Quasi-free bosons. In this case the Hamiltonian is given by H = l/^X^eA? 5 ? + 
l/2J2ijeA Vij x i x j> [ x i>Pj] — ifiSij.The harmonic coupling matrix V is symmetric, 
non negative definite, local on A and ||V|| = 0(1). The ground state is known to 
be a Gaussian state \V) with covariance matrix r(V) = V -1 / 2 © V 1 ! 2 and the gap 
with the first excited state is given by A = 2A m i„(V Al / 2 ) (A m j„(X) := minimum 
eigenvalue of X) [30J. The ground state fidelity can be obtained by performing 
a simple Gaussian integral and yields: T = \(V\V')\ ~ [det^- 1 + (r') _1 )] _1/4 ; 
setting r' := r(V + SV) = T + 5T and expanding at leading order in ST one 
finds for the fidelity susceptibility: \f ~ 1 1 r 1 t5T 1 1 ^ - This quantity can be bounded 
above as follows: X f < ^ / Kn m ( T )\\ 5T \\l = o(\\ST\\ 2 / A 2 ). Moreover, from the 
expression of T(V) above, one obtains ST = 1/2{V- 1 / 2 5V 8 V~ 3 / 2 SV) whence 
||5T|| 2 = lMWV-WdVh + ||^- 3 /W|| 2 ) = 0(A- 3 / 2 ||^|| 2 ). Finally, XE < X f = 
0(A- 5 rank(5T/)||5^|| 2 ) = 0(\dA\). In the last steps we exploited Eq. © and the fact 
that the boundary perturbation SV has rank 0(|cL4|) and O(l) norm. 



5. Area law corrections in gapless systems 

In this section we would like to illustrate how the entanglement susceptibility @ is 
able to capture the modifications to the strict area law that one expects in critical 
i.e., gapless systems. To this aim let us consider a d-dimensional tight-binding model 
on a hypercubic lattice. We partition the system into two half-spaces A and B of 
width La — 1 and Lb — 1 separated by a (d — l)-dimeiisional hyperplane and use 
the parametrization x = (xu,x±). To be consistent we fix the boundary conditions 
to be periodic along xii and open along x±. For simplicity we fix the size of all the 
parallel directions to L. The boundary term has the form: Hg = V + where V := 
Sx||G<9 c a!|[ rr° C x\\ x° an< ^ we can se ^ x< ± = ^ without loss of generality. The boundary 
operator V can then be expressed in terms of the Fourier basis which diagonalize Ha/b 
and turns out to be V = y/4/(L A L B ) Y,k v k ± , g± sinC^O sin (Qi-) c hl g± c k h k ± - Given 
the form of the one particle dispersion efc = 2 cos(fci) the energy denominator in 
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eq. ([2]) depends only on the perpendicular variables. The entanglement susceptibility 
takes the form 

8 ^ sin 2 (fcj_)sin 2 (g_L)n fc (l - n q ) 
L A L B , ^ {e k - e q y 

where the occupation numbers are given by n k = z9( — efc^ ) is the Heaviside step 
function) and we used the notation q = (k^,qj_). For large L we can introduce the 
function E(k±,q±) = J dk\\ n,k«,k±(l — nk»,qx) analogous to E(q) defined in [32]. 3 
measures the volume of feii such that (kn , k±) lies in the Fermi sea while (fey , q±) does 
not. For the specific tight-binding model at half filling S can be computed exactly in 
terms of a difference of volumes of the standard simplex. In general, for small 77 := q— k. 
one has 3 (k, k + 77) cx 7779(77 — o~k —Cq) where o- q //. = n jL A j B are the infrared cutoffs in 
momentum space. For large sizes one obtains \e = L d ~ x JJ^ dk J^ h+a dr) f {k,k + rf). 

Expanding f(k,k + r)) = /_i (k) jf] + YlnLo fn (k) T] n one realizes that the only 
divergent contribution is given by the integration of f-i(k)/n, while all other terms 
are convergent as (Jk/q 0. The integration over 77 gives trivially a log, and the 
diverging term turns out to be \e ~ —ctdL^ 1 ln(<7fc + a q ) where the constant is given 
by ad = f Q f-i (fc) dk [31]. If both La and Lb scale with L we obtain, for large L, 

XE-L^lnL. (8) 

This result is consistent with the logarithmic violations to the area law for fermionic 
systems discussed in [32] and [33]. We would like to stress here the two essential 
ingredients needed to obtain Eq. ([8]): i) linearity of 3 (fc, k + rj) for small 77, more 
precisely 3 (fc, A; + rj) oc 771? (77 — Uk — c<j); and ii) singularity of the form (q — k) for 
the energy denominators in Eq. ([7]). These features are believed to be general for 
free systems and rely on the existence of the Fermi surface and the linearity of the 
particle-hole excitations close to the Fermi energy. Similar assumptions have been 
used also in [32]. We thus expect Eq. ([7]) to be valid in more general cases with more 
complicated geometry and band structures, at least when the above assumptions are 
satisfied. 



6. Beyond the perturbative regime 

Low energy scaling behavior of a physical quantity is a property of the universality class 
(in the renormalization group sense) of the Hamiltonian. As long as the perturbation 
is not strong enough to induce a quantum phase transition scaling behavior is not 
expected to change. In other words, even though our argument for the entanglement 
area law relies on a perturbative expansion we expect it to hold true beyond the deep 
perturbative regime i.e., all the way up to A = 1 in Eq JT]). To elaborate further 
on this point let us for the moment come back to consider the ground state fidelity. 
The quantity log.F between the two ground states of the regions A and B when the 
boundary coupling Hg is switched on and off has been considered in [21] and named 
bi-partite logarithmic fidelity. The authors of [22J provide a heuristic argument, based 
on the standard quantum mechanics- classical statistical mechanics correspondence, to 
show that the bi-partite logarithmic fidelity fulfills an area law. The results discussed 
in this paper so far are tantamount to a rigorous proof of that claim to the second 
order in A. 



Entanglement susceptibility: Area laws and beyond 



7 



To investigate higher order contributions let us notice that for a gapped system 
with non-degenerate ground state one has T(X) = lim^oo N(f3)/N 1 ^ 2 (2(3), N(/3) :— 
(e~P H (^) where (• • •) denotes the quantum-mechanical average over |^ (0)) [21] . Us- 
ing the standard interaction picture formula e~^ H — e~P H °T s exp ^— A dsHg(s) S j , Hq(s) := 
e sH ° Hge~ sH ° one can prove [24^ 

(T s exp ( -A f dsH 9 {s) \ ) = exp ( ^(-A)"c„(/3) 




S„-l 



08)=/ d Sl --- ds n (H d ( Sl )---H d (s n )) c 



where the subscript c denotes connected (imaginary time) correlations functions. One 
can then write: 

log J-(A) = lim f>A)" (c„(/3) - \c n {2p)\ . (9) 

The presence of a spectral gap in H implies that connected averages scale as /? 
or more precisely c n ((3) = A n f3 + B n + 0(e^ l3A ) (A n ,B n are (3 independent) [24] ; 
whence log T(\) = i Y^=2 (~ ^)™ Bn- The term B 2 is easily seen to yield the fidelity 
susceptibility (B2 — ~Xf) and therefore, as we have seen, it is upper bounded by \dA\ 
in the gapped case. 

The key point is now that all these remarks concerning fidelity are directly relevant 
to Renyi entropy itself. Indeed it turns out that Renyi 2-entropy is a particular 
instance of logarithmic ground state fidelity [25]. To see this fact let us write the 
purity in the form Trp^ = Tr [5 , i 3 |4'o(A))(4'o(A)|® 2 ] where 613 is the swap operator 
between the first and third factor in (Ha £3> T~Lb)® 2 i-e., the state space of two copies 
of the system [23J. This trace can be viewed as the ground state fidelity between the 
ground state |* (A))® 2 of := H(X) ® 1 + 1 ® H(X) and 5i 3 |*o(A))® 2 i.e., the 
ground state of S13H ^Si3 [25\. Analog constructions extend to all a-entropies with 
integer a |25J. The fidelity formula © shows that one has a direct representation 

of the Renyi entropy in terms of connected (imaginary time) correlation functions 

(2) (2) 

of the boundary perturbation Vg :— Hg ' — Si^Hg 'S13 in a "doubled and twisted" 

theory. In particular all the entanglement susceptibilities -~ d n S2(0)/d\ n can 
be expressed by analogous B%" lst terms as defined above, associated with Vg [26] . 
Note that, for finite range interactions, Vg has O (\dA\) number of terms and in 
general represents a surface term. Hence extending to this doubled theory the heuristic 
statistical mechanics argument of [22] suggests that the area law should indeed hold 
for all the a- Renyi entropies S a (p) = — (a — l) _1 Trlog j o Q , a £ N) in the gapped case 
with clustering |34| . A deeper analysis of this argument may pave the way towards 
the rigorous understanding of how area laws for quantum entanglement and their 
violations arise in physics [35]. 



7. Conclusions 



Area laws for quantum entanglement are one of the most common features showing 
the strong atypicality of physical quantum states in the Hilbert space. In this paper we 
introduced a natural perturbative object, the entanglement susceptibility ([2]). Roughly 
speaking, this quantity measures the rate of entanglement generation (as measured by 
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the Renyi 2-entropy) when an infinitesimal boundary interaction term is switched on 
between two previously decoupled regions. Entanglement susceptibility while much 
simpler than entanglement entropy itself is able to describe its scaling behavior. 
This can be seen by bounding entanglement susceptibility with fidelity susceptibility 
[TU1 H5] • Both area law bounds for gapped systems and possible corrections for critical 
ones can be then obtained by elementary means. We illustrated the area law bounds 
for entanglement susceptibility with explicit calculations for the quasi-free systems 
and logarithmic corrections in the free fermionic critical case. 

The non-perturbative regime can also be explored by realizing that Renyi entropy 
is a particular instance of ground state fidelity in a doubled space |25| . This unification 
allows one to express entanglement by an exact series expansion of integrated 
connected correlations of a twisted boundary interaction. 

Let us conclude by mentioning a couple of natural goals for future research. On 
the experimental side, very recent work indicates that entanglement susceptibility 
may be directly measurable in bosonic optical lattices [33]. On the theoretical side, 
we would like to extend the results of this paper to excited energy eigenstates of local 
Hamiltonians [37] : Are they going to be more alike typical quantum states? 
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Appendix 

Proof of Let us expand the ground state of ([T]) in powers of A : 

|* (A)) = NJ2n=o Xn \y (n) ) d* (0) ) = I°a)|0b).) The perturbative eigenvector 
corrections |\I/( n )) are given by the elementary perturbation theory: = 
G Q (E ){H 9 \y( n -V) - J2kZo E[n ' k) \^ {k) ))i where £ (fc) =k-th order correction to the 
unperturbed eigenenergy E<® = E , G := J2 n >o( E o ~ £n) _X |*n)(#T»| = Q{E - 
Hq)~ 1 Q is the (projected) resolvent of Hq (Q ;= 1 — |>&o)(^o|0 The reduced density 
matrix is given by p A {\) = A/" 2 £^° =0 A V n) where pW ■- f^ =0 Tr B \^ {n ~ k) ){^ {k) \- 
Let us now consider the leading contributions n = 0,1; performing the partial trace 
over the region B one finds 

P {0) = \0a)(0a\, P (1) = ]T C(p A ,0)\0 A ){p A \ +h.c, 

PA>0 

P {2) = J2 \pa)(Qa\J2 C (p^Pb) C (<Ia,Pb). (.1) 

PA,qA PB 

in which C(pa,Pb) — {paPb\Hq\Qa®b) / {E PA + E PB — E ). Let us now compute 
the purity of p A : Trp 2 A = Af 4 Tr ((p^) 2 + \ 2 (p^) 2 + A W 2 '}) + 0(A 3 ). Here 
one has to notice that Tx{p^\p^} = 0, Ti{p^\p^} = 2J2 Pb> q\C{0 iPb )\ 2 , and 
Tr(pW) 2 = 2 ^2 pA>0 \C(pa, 0)| 2 . The normalization factor at this order is given by 
N- 2 = |||*o) + |* (1) )|| 2 = 1 + A 2 E Pa+Pb >o \C(pa,Pb)\ 2 - Using this explicit form and 
expanding again at leading order one obtains Trp A = 1 — 2 A 2 J2 PA PB >i \C(pa,Pb)\ 2 + 
0(A 3 ). Finally taking the negative of log one finds, the expression □ 
This derivation shows that the appearance of the fidelity susceptibility in 
our entanglement argument while remarkable is certainly not surprising. Indeed 
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an important ingredient in the derivation above is the the normalization factor 
1/Af 2 — || X)^Lo ^ n \^^) II 2 • Since in standard Rayleigh-Schrodinger perturbation 
theory one can always choose (v]/( ) |vj/( fc )) = for fc > 0, it follows that M 2 is 
nothing but the ground state fidelity: T 2 = |(* (0)|*o(A))| 2 = | (* (0) |*o(A)> | 2 = 
A^ 2 IE^o(*o|* (fe) )| 2 -AA 2 . 



